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^H ! Abstract 

^*^ . A numerical and analytical study of the relaxation to equilibrium of both the Fcrmi- 

G ' Pasta-Ulam (FPU) a-model and the integrable Toda model, when the fundamental 

mode is initially excited, is reported. We show that the dynamics of both systems is 
almost identical on the short term, when the energies of the initially unexcited modes 
grow in geometric progression with time, through a secular avalanche process. At the 
CN ' end of this first stage of the dynamics the time-averaged modal energy spectrum of 

the Toda system stabilizes to its final profile, well described, at low energy, by the 
spectrum of a g-breather. The Toda equilibrium state is clearly shown to describe well 
i^Q . the long-living quasi-state of the FPU system. On the long term, the modal energy 

CN ' spectrum of the FPU system slowly detaches from the Toda one by a diffusive-like 

r«^^ , rising of the tail modes, and eventually reaches the equilibrium flat shape. We find 

l^ ' a simple law describing the growth of tail modes, which enables us to estimate the 

time-scale to equipartition of the FPU system, even when, at small energies, it becomes 
unobservable. 
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1 Introduction 

The present work deals with the so-called Fernii-Pasta-Ulam (FPU) problem, which con- 
sists in understanding and characterizing the process of dynamical relaxation to the micro- 
canonical equilibrium of non integrable Hamiltonian systems with many degrees of freedom, 
when initial conditions far from equilibrium are chosen. This issue, under investigation for 
more than 50 years, is relevant to equilibrium and out-of-equilibrium statistical mechanics, 
as well as to many problems in condensed matter physics. 

Following the original work of FPU [1], one of the most studied systems has been, and 
presently is, the so-called a-model, defined by the Hamiltonian 
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where q is a parameter and fixed ends conditions are chosen: qo = Qn = Po = PN = 0. 
Such a Hamiltonian model with A^ — 1 degrees of freedom can be regarded either as a 
space-discretized version of a nonlinear string (as FPU did), or as the expansion around 
the equilibrium configuration of a one-dimensional system of pairwise interacting identical 
particles, and is commonly referred to as a nonlinear oscillator chain. 

FPU numerically integrated the equations of motion of system ([1]), given by (jn = 
dHa/dpn, Pn = —dHa/dqn, and they chose initial conditions of the form 

g„(0)=Asin(^) , PniO) = 0, (2) 

A being a parameter, thus initially exciting the fundamental (longest wavelength) Fourier 
mode of the chain. FPU considered system sizes N = 16, 32, 64. The surprising and for a 
long time debated result they got, namely the absence of any complete relaxation of the 
system to a state of equilibrium compatible with the laws of statistical mechanics (at least 
up to the observation times then available) was named, after them, the FPU paradox. As 
first conjectured in Refs. [21 [3] and then clearly shown in Ref. [4], the paradox is due 
to a separation of the relevant time-scales of the problem, which in turn depends on the 
value of the total energy E and on the size A^ of the system. More precisely, it has been 
understood that energy is first transferred to shorter wavelength modes, which, within a 
first characteristic time-scale T'^^{E,N) gives rise to a quasi-state (as named by FPU, later 
also referred to as "metastable state" [21 [3] or "natural packet" [1] in the literature) . This 
state is characterized by a modal energy spectrum displaying an exponentially decreasing 
tail. On a second, larger time-scale T'^''{E,N), the modal energy spectrum approaches an 
almost flat shape, i.e. the system does reach the equilibrium state characterized by modal 
energy equipartition, as predicted by the laws of equilibrium statistical mechanics. The FPU 
paradox appears then when, depending on E and N, a strong inequality T'^^ <^ T^'^ holds, 
such that the approach to equilibrium becomes unobservable on the available computation 
times. 

Up to now, the short term dynamics of the system has been thoroughly investigated, 
both from a numerical and an analytical point of view, also for initial conditions more 
general than ([2]), and both the time-scale T''^ and the slope of the exponential tail of the 



energy spectrum have been more or less quantified. In particular, the localization of the 
energy in Fourier space, i.e. the existence of the quasi-state, has been explained essentially in 
two ways. One approach consists in looking for a particular stable solution of the equations 
of motion that is close to the actual solution of the initial value problem and displays 
a localized energy spectrum. One way to get such a solution consists in looking for the 
Lyapunov continuation of the initially excited mode, or g-breather [3 [U O [HI [9] . The other 
approach consists instead in looking for a suitable integrable system whose dynamics is 
close, on a certain time-scale, to that of the actual system, and leads to the Korteweg-de 
Vries equation and its truncations in mode space [9l [lOl HH EZl ESI HU US] (Ref. [9] contains 
a first tentative comparison between the two approaches). On the other hand, for the a- 
model ([T]), there exists an approximating integrable system given a priori, namely the Toda 
model (see Hamiltonian ([3|) below), and a comparative study of the two models has been 
performed e.g. in Refs. [H [HI [HI H [20] . 

Notwithstanding the results just quoted, both a systematic study of the dynamics of 
problem ([I])-(l2|) on short times and a deeper understanding of the possible existing links 
between the aforementioned theoretical approaches are still lacking. Such a gap is partially 
covered in the present paper. We explain the detailed resonance mechanisms ruling the 
energy cascade on short term together with a simple characterization of the quasi-state. 
Moreover, we compare, for the same initial conditions of the form ([2]), the dynamics of the 
FPU chain with that of the integrable Toda modelpT] defined by the Hamiltonian 
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with fixed ends at n = and n = N . The integr ability of the latter Hamiltonian system 
was first suggested in Ref. [22] and then proved in Refs. [23] and [23]. Expanding the 
exponential in ([3|) and taking into account the boundary conditions, one can write 

N-l , _ .^ 

n=0 r->4 

i.e. the FPU Hamiltonian can be regarded as a perturbation of the Toda one. The call for 
such a comparative study for the evolution towards an equilibrium comes from the necessity 
to distinguish those phenomena that are of integrable - i.e. non chaotic - nature, from those 
that are instead signatures of non integrability - i.e. chaos. In the present paper we show 
that the FPU quasi-state coincides well with the Toda equilibrium state (see below. Sections 
II and HI; see also Ref. [16]), and is therefore a manifestation of closeness to integrability. 
Chaos becomes relevant in the a-model only after the saturation to the quasi-state. We 
also show that the quasi-state, i.e. the Toda equilibrium state, is well described, at low 
energy, by a time-periodic orbit or one-dimensional torus, namely the g-breather associated 
to the first mode. 

For what concerns the second time-scale T^^, its numerical estimate has been given 
e.g. in Ref. [18] for the a-model, and in Refs. [25l [26l [27] for the so-called /3-model; the 
two models are treated together and compared in the more recent Ref. [28] , where accurate 



numerical estimates of T^'^ are given. Yet, a detailed description of how the tail of the modal 
energy spectrum raises from zero to the equipartition level has not been provided so far. 
In the present paper, we perform a fit of the numerical data relative to the modal energy 
spectrum of the a-model during the second part of its evolution, which yields a simple 
law describing the motion of the tail. This allows us to measure the second time-scale T'^'^ 
on times which can be much shorter than T^'^ (notice that a direct measurement of the 
equipartition time-scale T^'^ becomes impossible when it exceeds the available computation 
times). 

We find that, at low energy (when the quasi-state is a one-dimensional Toda torus) the 
estimated time-scale T^^ undergoes a transition from power law to stretched exponential 
of the inverse of the energy. As a consequence, the quasi-state of the a-model becomes 
effectively stable: in practice the approach to equipartition becomes unobservable, which 
explains the FPU paradox. 

Let us stress that the existing literature on the FPU problem is huge, encompassing 
more than fifty years of research on an unsolved problem. We thus defer the reader to some 
of the existing high-quality reviews on the subiect[29 l[30l[3T] . 

The paper is organized as follows. In Section II, after introducing some notations and 
a few relevant quantities, we state our main results on the evolution towards equilibrium 
of systems ([T|) and ^ with initial conditions ^. In Section III, the phenomenology of the 
problem is illustrated by the numerical results and the fits with the analytical predictions. 
In Section IV, the theory leading to the analytical estimates is reported. Finally, Section V 
is devoted to concluding remarks. 

2 Main results 

We integrated the equations of motion associated to the Hamiltonians ([1]) and ([Sj) by means 
of a symplectic algorithm, namely the Yoshida kinetic-potential splitting algorithm of fourth 
order [321 [33], with a time-step kept fixed to 0.05 all over the explored energy range. [M] 
The size N of the system considered in the present paper is kept fixed to A^ = 32, and no 
systematic exploration of size dependent phenomena is performed. The chosen value of the 
nonlinear ity parameter is a = 0.33. [35] 

We follow the evolution of the harmonic energies E^ of the Fourier modes of the system, 
also coined modal energies. These quantities are defined as 

^^(,) ^ ^(t)+|M(!) , (, = i,...,Ar_i), (5) 

where {Qk,Pk) = a/2/A J2n=i {Qn,Pn)sin{iTkn/N) are the Fourier coordinates and 

ujk = 2sml—j (6) 

is the dispersion relation of the linearized problem (a = 0). For the FPU model ([1]), the 
total energy Ea = Ha{q{0),p{0)) in terms of the initial condition ([2]) reads 



whereas, for the same initial condition, the total energy Et = H'T{q{0),p{0)) of the Toda 
model ([3]) reads [36] 

where loix) is the first modified Bessel function of order zero |37j . All the comparative runs 
were performed at the same value of the total energy E for the two models, which required 
to consistently choose different initial amplitudes A and A' in the two cases: E = Ea{A) = 
Et{A'). [38] 

For a = the modal energies E^ do not evolve in time both for the FPU and the Toda 
model, since they are constants of motion: Ei = E and Es>2 = 0, for all times. When 
a > the modal energies are no longer constants of motion for both systems, and evolve in 
time. Now, for the nonintegrable FPU system, one expects that each time-averaged modal 
energy 

Ek{t) = - / Ek{s) ds (9) 



converges, as t — )• oo, to its final expectation value [39]. which is approximately equal to the 
specific energy E/N of the system. 

On the other hand, for the integrable Toda model the modal energies are quasi-periodic 
functions of time (with a maximum number A^ — 1 of rationally independent frequencies), 
which ensures that their time-averages @ converge to some limit as t — )■ oo. Such a limit 
defines what we refer to as the equilibrium state of the Toda system. In the sequel we use 
and plot the normalized modal energies 

e.(t) = ^ (10) 

and their time-averages efc(t) = Ek{t)/E, and we refer to the plots of e^ vs. k and of e^ 
vs. k as the instantaneous and averaged modal energy spectrum (m.e.s.), respectively. The 
initial values of the normalized modal energies are then ei(0) = 1 and ei(0) ~ 1 — a^E/N , 
for the FPU-Q and Toda model respectively, whereas 6^(0) = for k = 2, . . . ,N — \. The 
expected equipartition value of the e^'s is approximately 1/A^. 

We report in Fig. [1] nine snapshots of the FPU and Toda averaged m.e.s. at times 
t = 0, 10, . . . , 10^, all referring to the same run at E = 1. We note that up to time t = 10"^ 
the two spectra almost perfectly superpose, whereas beginning with t = 10^ the tail of the 
FPU m.e.s. starts to rise while the whole Toda m.e.s. no longer evolves. Finally, at t = 10*^ 
the FPU m.e.s. displays an almost perfect flat shape, a signature of statistical mechanical 
equilibrium. 

2.1 Stage I: secular avalanche 

On a very short time-scale, the initially excited first mode {k = 1) looses its energy, while 
the energies of all the other modes increase with a mode-dependent power of time. For both 
the FPU and the Toda models the time-evolution of the instantaneous m.e.s. is described 
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Figure 1: e^ vs. k (log-linear scale) for FPU-q ([T]) [orange triangles] and Toda ([3]) [green 
squares] models, ii^ = 1, A^ = 32, a = 0.33. From the upper left corner to the lower right 
one t = 0, 10, . . . , 10^. Note the change of vertical scale from t = 10^ on (last six panels). 



by the following law: 
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where 
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is a small parameter and the time-independent coefficients Ck are recursively computed by 
solving the nonlinear equations 

fc-i 
ci = 1 ; Ck = ^r^^^Ck-gCg {2<k<N-l) . (14) 

q=l 

These results follow from the smallness of the difference ujk — kuji, and essentially consist 
of a resonant transfer of energy, where the secular (proportional to t^) growth of the modal 
energy 62 in turn implies the law e^ oc t^''^"^) (A; > 2). We call this stage of the cascade 
the secular avalanche stage (see Section [3. ip . The relevance of the acoustic resonance as a 
driving mechanism of the cascade in the FPU problem was first stressed in Ref. [30] . 

2.2 Toda equilibrium and quasi-state 

Until a characteristic, energy-dependent time-scale 

N 
T"' « — (15) 

the decrease of ei and the resonant growth of efc>2 due to the secular avalanche effectively 
slow down in the FPU case and actually stop in the Toda case. The time-averaged m.e.s. 
of both the FPU and Toda models converge to an exponentially localized profile, which, for 
energy values up to about 10~^, is well described by the expression 

with ei determined by the normalization condition 

N-l 

fc=i 

The m.e.s. profile (J16p is that of the (7-breather associated to mode k = 1, i.e. the Lyapunov 
continuation of the first mode to the nonlinear system [5]. The single energy-dependent 
frequency of this periodic orbit, or one-dimensional torus, is 

f^i(^) ^u:, + M^ei /^^^^ + . . . (18) 

!i[ZU)i — U)2j 

Though the agreement between formula (I16p and the numerical m.e.s. becomes the less 
accurate the higher is the energy, we numerically observe that the m.e.s. of both models 
keep on saturating to (almost) the same exponentially localized profile for all the explored 
energy values up to £" = 1. Such an agreement suggests a simple interpretation to the 
nature of the quasi-state guessed by FPU: the quasi-state of the a-model coincides with the 
final equilibrium state of the Toda model. It must be stressed that strong indications to 
such a conclusion can be found already in Refs. [161 120] , 

The analytical predictions ()lip - (|12p and (J16p . were obtained by computing a particular 
resonant normal form Hamiltonian of the system to leading order. The same normal form 
is the starting point of heuristic considerations concerning the observed stability of the 
q-hreathev m.e.s. (I16p on long times. 



2.3 Stage II: diffusion of tail modes 

After saturation, the FPU ni.e.s. slowly continues to evolve and detaches from the reference, 
stationary Toda ni.e.s. . This behavior becomes actually unobservable below some energy 
threshold Ec {Ec — 0.1 for N = 32), at least up to the available observation times (lO*^ to 
10^). Above this threshold, there is an evident raising of the m.e.s. tail. We consider the 
evolution of the normalized total energy residing in the last third of the modes, namely 

31 

7?(t) = Y. ^'^(*) ' (19) 

fc=22 

The choice of the last third of modes in defining the tail of the m.e.s. is of empirical nature, 
and agrees, when the quasi-state is well described by the g-breather, with the definition of 
tail given in Ref. [8] and based on resonance arguments. We find that, after a transient 
time approximately coinciding with the time-scale T'^^ , the quantity rj{t) starts to increase 
with time as a power law 

m ~ Dt^ , (20) 

within an energy depending time-window. Both the exponent 7 and the coefficient D depend 
on the energy. The diffusion exponent '~f{E) roughly displays a step-like behavior: 'y{E) ~ 
for E <^ Ec and 7(-E) ~ 1 as S S> Ec- The energy dependence between the two limits 
is not smooth: the exponent ^{E) strongly fluctuates around Ec, the latter value being 
more properly identified as the center of a transition interval rather than an actual sharp 
threshold. This is possibly a signature of anomalous diffusion processes, characterizing 
the dynamics of tail modes (recall that e^ = Ei^/E depends quadratically on the modal 
coordinates Pk and Qk)- Depending on the value of the energy E one can observe either 
sub-diffusion (7 < 1) or super-diffusion (7 > 1). An analytical derivation of the law ()20p is 
currently missing. 

We then extrapolated the time-scale to equipartition T^'', defined as the time necessary 
for its time average rj to reach its approximate equipartition value 1/3: 

T"" = {3Dy^^^ . (21) 

Both D and 7 can be therefore obtained on relatively short times. These values lead to lower 
bounds for the actual equipartition time, which for example can reach the unconceivable 
value of 10^°°° at £; = 0.01. We find that T^^ oc S"" at higher energies (from E = 0.5 up 
to 2), with a close to 3, whereas T'^'^ ~ exp(c/£^^) for low energies (in the range 0.01 to 0.1), 
with 6 ~ 0.7. 

In this second stage of the dynamics we also measured the erosion of one of the first 
integrals of the Toda model independent of the Hamiltonian, computing it along the FPU 
dynamics. It turns out that, above the energy threshold Ec, when the diffusive-like raising 
of the tail modes becomes observable, such a quantity undergoes a drift starting at times 
of the order of T'^^. This behavior indicates that the relaxation to equilibrium of the FPU 
system is characterized by the drift of the Toda integrals, which become quasi-invariants of 
the FPU system at low energy. 



3 Phenomenology and Numerics 

We now come to a detailed description of the phenomenology of the problem, presenting 
our numerical results and their comparison with the theoretical predictions. 

3.1 Stage I: secular avalanche 

The secular avalanche stage is shown in Fig. [21 where the log-log scale plot of efc(f) vs. 
t at E = 0.01 is reported, both for the Toda and the FPU-a system, together with the 
lines corresponding to the theoretical law (|lip - (jl2p . The values of the mode number are 
restricted to fe = 1, . . . , 10 just to improve the readability of some details in the figure. One 
can note the almost perfect superposition, mode by mode, of the Toda and FPU numerical 
modal energies, which are both well described, on the very short term, by the law (|11|) -(|12|). 
The same good agreement of the latter theoretical law with the numerical data is found all 
over the energy range 10~^ to 1, though further figures are not reported. 

The law efc(i) oc f^C^'"^) for low modes is due to the smallness of the difference kuji — ujf^, 
i.e. to their almost complete resonance. Indeed, mode A; = 2 is forced in resonance by 
the second harmonic of mode A; = 1, so that its amplitude grows linearly and its energy 
grows quadratically with time: 62 oc t^. Mode /c = 3 is in turn forced in resonance by the 
product of amplitudes of modes k = 1 and k = 2, which has a factor t in front, so that 
ek oc (t^)^ = t^, and so on. In this way a resonant cascade of energy transfer sets in, which 
causes the growth of the initially unexcited modes one after the other. Notice that for high 
modes, the difference w^ — kcoi is not small, but their amplitude is initially zero, and this 
(as shown in Section IV) ensures the validity, also for them, of the approximation leading 
to the secular avalanche law (jlip -(|12p. 

Let us denote by T^'^{E,N) the mode-dependent characteristic time of validity of the 
secular avalanche law (jlip - ()12p . For the first mode this time can be estimated by simply 
requiring that the leading order estimate (jlip yields ei > 0. Setting ei = in (llip one 
gets Tf"- = 2\/^/{fi^LOiuj2) — 2N/{tt^). Moreover, in the regime of weak mode-coupling 
explored here, almost all the energy is retained by the harmonic part of the Hamiltonian, so 

that the sum X]fc=i ^k is close to its initial value (one) for all times, which in turn implies 

1 

efc < 1 for all /c's. Setting e^ = 1 in ([12]) one gets T^'^ = c^ ^~^ j [i for fc > 2. We finally notice 
that T^ ~ T'^" = 0.2 N j [i for large values of k. By exploiting such asymptotic behavior, 
the law ([12]) reads, for A; > 1, efc(t) = e"^^''"^)'^^*), where a{t) = -ln(t/r^"). A logarithmic 
time-dependence of the slope —cr{t) of the m.e.s on the short term was first predicted, with 
methods based on continuum approximations, in Refs. [21 [JT]. It is interesting to notice that 
the time-scale ()15p . derived below (in Section [4. 5p from the point of view of mode dynamics, 
was first determined as the break-down time of the solution of the nonlinear wave-equation 
approximating the FPU dynamics in the continuum, zero dispersion limit. [42] 143 j 

3.2 Toda equilibrium and quasi-state 

At the end of the secular avalanche, the time-averaged m.e.s. (e^ vs. k) of both the FPU-a 
and the Toda models saturate to one and the same exponentially localized profile - the FPU 
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Figure 2: efc(t) vs. t (log-log scale) for A; = 1, . . . , 10 at -E = 0.01; black and orange curves 
refer respectively to the Toda and the FPU numerically computed energies, while blue lines 
correspond to the theoretical law (llip -()12p. 

quasi-state. For larger times the Toda time-averaged m.e.s. no longer evolves, whereas the 
FPU one keeps on evolving, the more slowly the lower is the value of the energy E. In order 
to illustrate this, in Fig. [3] we report both the time-averaged Toda and FPU m.e.s. and 
the FPU instantaneous one (log-linear scale). The nine panels refer to three values of the 
energy £^ = 0. 01, 0. 1,1 at times t = 10^, 10^, 10*. The value of the energy E is fixed in each 
row, whereas that of the time is fixed in each column. At E = 0.01 (first row) the Toda and 
FPU spectra are almost perfectly superposed at any time, the only observable difference 
being the small fluctuations e^ — e;^. of the FPU energies of the last (highest) few modes. At 
E = 0.1 (second row) the FPU spectral fluctuations 6^ — 6^ become observable for any k, and 
increase with time. Also at t = 10^ the tail of the FPU m.e.s. starts to rise detaching from 
that of the Toda. Finally, at E = 1 (third row) (see also Fig. [T]) one observes, in the FPU 
m.e.s., even larger fluctuations, the raising of the tail and the relaxation to equipartition. 
Fig. [3] also shows that at low energy, when the FPU model takes a long time to relax to 
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Figure 3: ek{t) vs. k (log-linear scale) for the Toda (green squares) and the FPU-a (orange 
triangles) systems, and efc(i) vs. k for the FPU-a model (blue stars). E = 0.01,0.1, 1 from 
top to bottom, t = 10^, 10^, 10^ from left to right, i.e. the energy is constant along the rows 
and the time is constant along the columns. Notice the change of vertical scale from top to 
bottom. 



equilibrium, its motion stays close to the invariant torus of the Toda model corresponding 
to (almost) the same initial condition. The slow variation along the FPU actual trajectory 
of the Toda invariants is at the basis of the FPU paradox. This was first conjectured in Ref. 
|16| . where a phenomenological picture of the Toda dynamics based on inverse scattering 
is presented. In particular, in the latter reference, it is shown that at small energy, the 
Toda dynamics corresponding to the initial excitation of a single mode takes place on a 
one-dimensional torus, and that such a periodic orbit must be the Lyapunov continuation 
to the nonlinear problem of the linear normal mode solution. 

With this in mind, in Fig. Uthe time-averaged m.e.s. (in log-linear scale) of the Toda 
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Figure 4: Time-averaged Toda m.e.s. at three different times (10"^, 10^, 10^), E 
asterisks correspond to the theoretical prediction (J16p . 
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model at E = 0.01 is reported for three different values of the time t = W^, 10^, 10*, using 
different symbols. Notice that the three numerical spectra at different times almost perfectly 
superpose, which shows that the time-averaged m.e.s. we are looking at is stationary. The 
theoretical m.e.s. (fT6]l (with ei computed from (fT7|l ) is also plotted in Fig. [H (asterisks). 
The agreement with the numerical data is very good and holds for energies lower than 
E = 0.01. On the other hand, in Ref. [16] it is also shown that the dimension of the 
Toda torus corresponding to the a single mode excitation grows with the energy, so that 
the (Toda) motion becomes quasi periodic. This explains why the g-breather m.e.s. (|16p 
fails to describe well, at larger energy values, the numerical spectra of the Toda and FPU 
systems, although the latter spectra go on to almost superpose on short times. The analytic 
description of the Toda or FPU m.e.s. corresponding to the full quasi-periodic case is beyond 
the scope of the present paper and is left to future investigation (see however Refs. |44^ I45j 
for interesting results in this direction). 
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3.3 Stage II: diffusion of tail modes 



10^ 



"T 1 1 I I I 1 1 1 I I I 1 1 1 I I I 1 1 1 I I I 1 1 1 I I I 1 1 1 I I I 1 1 l—TT 



A 'i ^/^v" ' ^y"" ' Y^ ' * ' ^" 



E=0.1 - 




Figure 5: Bottom to top: rj{t) vs. t (log-log scale) for E = 0.01 (blue) 0.1 (magenta) and 
1 (green), together with the corresponding time-averages rj(t) vs. t (black curves) and the 
fitting lines Df^ vs. t (dashed lines). 

In order to describe the approach to equipartition in the FPU system, in Fig. [5] we 
plot (in log-log scale) the instantaneous normalized tail energy r/(i) = Ylk=22 ^k{t) and the 
corresponding time-averaged values r/(t) vs. t, for three values of the energy E = 0.01, 0.1, 1 
(bottom to top). Note that the instantaneous values evolve by fluctuating around the 
corresponding time-averages. We observe a growth of both quantities, well fitted by a 
power-law rj^t) ~ Dt"' within some time-interval [to,ti]- It follows the linear behavior 
logT]{t) ~ logZ? -|- 7logt in log-log scale, as reported in Fig. \5\ Notice that when E = 1, r] 
and r] asymptotically reach the equipartition value 1/3 and do not grow further, as expected. 

By a systematic and detailed numerical investigation we first determined the time inter- 
vals [tQ{E),ti{E)] and then computed jiE) and D[E) by a least squares fit (200 values of 
E were considered, in the energy interval 0.01 to 2, with an energy step iS.E = 0.01). The 
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Figure 6: 7 vs. E corresponding to tj. 

errors on both quantities turn out to be small (a few percents of the value) except for very 
low energy values, close to 10~^, where large errors affect the almost vanishing exponent 7. 
The function ■~f{E) is plotted in Fig. [6l Despite strong fluctuations, we find that 7 changes 
from a value very close to zero at energies lower than 0.1 to values oscillating around 7 = 1 
for energies larger than 0.1. Moreover, the absolute value of the fluctuations decreases with 
increasing energy. 

We conjecture that the fluctuations in the exponent 7 are intrinsic to the dynamics, and 
are due to the complicated phase space structure of the FPU system. To show that, we plot 
r]{t) in Fig. [7] for two values of the energy which differ by AE = 10~^. Indeed, the growth 
rates differ substantially: the values of the exponent 7 corresponding to these two energies 
are 7(0.14) = 0.03 and 7(0.13) = 1.1. 
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Figure 7: r/(t) vs. t for two values of the energy differing by AE = 10 ^: the upper and 
lower curves correspond to the energy values E = 0.13 and E = 0.14, respectively. 

3.3.1 Estimate of the time to equipartition 

The observed power-law growth rj ~ Dt"' allows us to define the extrapolated equipartition 
time such that DP equals 1/3, i.e. the equipartition value of tJ, which yields T'^'^ = [3D)~^'^ . 
In Fig. [8]we plot log(r'^'') vs. E. We observe a sharp crossover around E = 0.1, below which 
j'eq increases dramatically, up to values T^'^ ~ lO^*^*^*^. Such large times strongly suggest 
that below the reference energy threshold Ec ~ 0.1, the quasi-state of the FPU system (i.e. 
the equilibrium state of the Toda model) might become stable in the sense of Nekhoroshev, 
i.e. over exponentially long times of the order 10"^' ^ , with suitable constants c and b. We 
get a rough estimate for the stretching exponent b ~ 0.7 it 0.2. Existing numerical results 
in the literature [IHl [261 EH EH] suggest one that for larger values of the energy the time to 
equipartition grows inverse proportional to a power of the energy E. In the inset of Fig. [8] 
we plot log(T^'^) vs. \og{E) in the energy interval 0.5 to 2. We indeed find some evidence for 
a power law, with exponent a ~ 2.9. This number is close to the value a = 3 first reported 
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Figure 8: Estimated log(r^'') vs. E from average 77 data. Inset: log(T^'') vs. log{E). The 
line drawn has slope —2.9. 



in Ref. Il8j for N = 32, and there obtained with other methods. 



3.3.2 Evolution of a Toda integral 

We have studied the evolution of one of the additional constants of motion of the Toda 
model. The integral J{q,p) we considered is [23] 



iV-l 



^-E 



n=0 



Pn , {Ql + Qn+iqn + ql+l)fn ^ {fn+1 + fn + fn-l)fn 



+ 



+ 



a 



^o? 



(22) 



where /„ = e^"('^"+i~'5'"V(2a) {n = 0, ...,iV - 1) are the Flaschka variables [21], with 
boundary conditions /_i = /o, Jn = Jn-i, inherited by the restriction of the 2A^-periodic 
case to that of fixed ends. Since J is an integral of the Toda model, J = {J, -ffr} = along 
the Toda flow. Along the FPU flow the quantity J = {J, Ha} = {J, Ha — Ht} is expected 
to be small in the perturbative regime where the two models can be considered to be close 
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Figure 9: |J^°'^"(i) 
and £■ = 1 (green). 



jioda(o)| vs. t (log-log scale) at E = 0.01 (blue), £ = 0.1 (orange) 



to each other. Let us denote by J (t) and by J (t) the evolution of the quantity (I22p 
along the numerical FPU and Toda trajectories, respectively. 

In Fig.[9]the absolute error | J-^°'^"(t) — J-'"°''''(0)| vs. t is reported, for the three different 
values of the energy E = 0.01, 0.1, 1 (bottom to top), in log-log scale. After a transient time, 
the error stabilizes quickly at a practically constant value ranging from 10~^ at E = 0.01 
to 10~^ at E = 1. Thus the numerical algorithm reproduces the integrable Toda dynamics 
with very good accuracy. 

In Fig. [To] we plot the evolution of \J^^^{t) — J^°'^'^{t)\, for the same three values of 
the energy E = 0.01,0.1,1 (bottom to top), in log-log scale. For E = 0.01 or 0.1 this 
difference stays almost constant (up to the final observation time 10*), while for £^ = 1 it 
starts to grow at about t = 10^ to eventually saturate at a constant value. Therefore at 
low energies the FPU system can be regarded as a perturbation of the Toda one, the exact 
integrals of the latter being adiabatic invariants, or quasi- integrals, of the former. Note that 
at E = 0.01 and 0.1 no effective tendency to equipartition was observed up to t = 10*. On 
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Figure 10: | J^^^(t) - J^°'^'^(i)| vs. t (log-log scale) ai E = 0.01 (blue), E = {).l (orange) 
and E = 1 (green). 

the other hand, for E = 1 the tail of the FPU m.e.s. starts to effectively raise at t = 10^, 
when the drift of the Toda integral starts up. 



4 Theory 

Our theoretical approach is based on Hamiltonian perturbation theory. We construct a 
resonant normal form Hamiltonian of the FPU and Toda systems. The laws (Ilip -(jl2p are 
then derived by solving the normal form equations, to leading order, for short times. The 
law (|16p is instead obtained as the m.e.s. of a time-periodic solution of the normal form 
equations. The theoretical explanation of the diffusion law l\20\\ goes beyond our present 
capabilities. 
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4.1 Complex modal variables 

Let us consider the Haniiltonian ([T]) of the FPU a-model. We first perform a change of 
variables passing from the real space coordinates {q,p) S M^^ ^^^ to the complex modal 
variables u € C^~^. The complex modal variables u^ are defined in terms of the real ones 
by 

Uk = 7== , \Mt)\ =—E^=^k{t), (23) 

whose utility is self-explanatory. The FPU Hamiltonian ([1]), when expressed in the new 
complex variables Uk, reads 

7V-1 N-l 3 

fc=l fci,fc2,fe3 = l i = l 



W2 -Hi 

where a superscript asterisk denotes complex conjugation, and the mode-coupling coefficient 

^ki,k2,k:i — "ki+k2,k3 + "k2+k3,ki + "k:j+ki,k2 ~ ^ki+k2+k3,2N (25) 

has been introduced {Sn,m is the usual Kronecker symbol). We also recall that 

M = a/l (26) 



denotes a small parameter in the theory. In ()24p we have explicitly pointed out the quadratic 
part 7^2 and the cubic part Ti^ of the Hamiltonian. The Hamilton equations of motion 
Uk = {uk,^^} = —iuJkdTi/dul associated to the Hamiltonian Ti explicitly read 
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p,q=l 















ilk = -i^k ■"fc + 4 2^ ^k,p,q{up + Up){uq -h n*) . (27) 

p,q=l 

Finally, the initial condition ([2|) expressed in the new variables reads 

uk{0) = 6k,i . (28) 

4.2 Averaging 

The initial condition (j28p . at low energy, is expected to excite an oscillation of the first mode 
at a frequency f^i close to the unperturbed one, the difference ill — uji vanishing as /i — t- 0. 
Indeed, in the latter limit, equations ([27|) become linear and are solved by Uk = e~*'^^*(5fc i. 
Now, by substituting the leading order ansatz Uk = e~'^^^^6k,i in the right hand side of 
equations (p7|) . and making use of ([25]) . one realizes by a direct inspection that the mode 
A: = 2 is forced, among others, by a term oscillating at a frequency 2rii ~ 2uji ~ uj2- 
Thus, mode k = 2 responds close to resonance and its amplitude grows, which requires to 
modify the starting ansatz to Uk = Ae~^^^6k,i + Be~ '^^6k^2- Re-inserting the latter in 
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the right hand side of (j27p one finds that mode A; = 3 is almost resonantly pumped at a 
frequency SJli ~ Sui ~ oj^, and so on. Such a process leads to an almost resonant transfer 
of energy from the fundamental mode to higher frequency modes, with mode k oscillating 
at the frequency kQi. For this reason, we perform the time-dependent change of variables 
{u,u*) I-)- {v,v*) 

{u,u*) = <^'''\v,v*) : ufc = e-^^^i*t;fc (29) 

for k = 1, . . . ,N — 1. The above defined map $^1* is actually the 27r/Oi-periodic flow of 
the Hamiltonian 



hn,{u,u*) = V ( — -) 



uk\. (30) 



Notice also that all the functions of the modal energies e^ = litfcP, such as 7^2 and hfi-^, 
are invariant under the flow $^1*, i.e. Poisson-commute with /i^^. The new variables v 
evolve according to the Hamilton equations v^ = {vk,A} = —icokdA/dv'^ associated to the 
explicitly time-dependent, 27r/r2i-periodic Hamiltonian A: 

A{v,v*,nit) = n2{v,v*) - hnAv,v*) + n3{<^'''\v,v*)) . (31) 

According to the usual averaging theory |461 I47j. we replace the explicitly time-dependent 
Hamiltonian A with its time average A. This amounts to replacing ^3 o $^1*^ on the right 
hand side of (I3ip . with its time average 

n3iv, ^*) = ^J "' ^3 {^""''{v, V*)) ds . (32) 

We then obtain an averaged system, whose solution stays close to that of the original 
system over times inverse proportional to the small parameter. The averaged Hamiltonian 
A = 'H2 — hn^ + Ti-s of the system then reads 

and its associated Hamilton equations Vk = {vk,A} = —ico^dA/dvl, read 

(fc-l N-k-l \ 

y^ Vk-qVq + ^ '2Vk+qV*q , (34) 

9=1 9=1 / 

where, if /c = 1 or A; = A^ — 1, the first or the second sum on the right hand side are 
respectively absent. Such equations must be solved, for A; = 1, . . . , A — 1, with the initial 
condition 

Vk{0) = 5k,i . (35) 

Note that since the averaging (j32p is performed along the flow of the Hamiltonian h^-^ , 
defined in (pOJ) . then hfi-^ is a constant of motion of the averaged system: {/i^j,^} = 0. 
Thus also the quantity 

h.Av,v*) = j:(—)\vk\' (36) 

k=i ^ ^'' ^ 
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is a constant of motion of the averaged system. 



4.3 Solution of the averaged system 

The averaged equations (|34p contain the small parameter /j,. We set 



(37) 



the dependence on /i being left in the nonlinear frequency shift 5oJi = Qi — oji. We require 
that 6uji —7- as /i —7- 0, so that in absence of nonlinearity the initially excited mode will 
oscillate with its unperturbed frequency uji. It follows [38] 



6uJi{n) = Oi(;U) -UJi= ^/i^JzyO) 

The initial condition ffc(O) = 6k,i implies 

4')(0) = 6k,i6j,o ■ 



(38) 
(39) 

(40) 



Upon substitution of ()38p - ()39p in (j34p and balancing order by order one gets a nonlinear 
system of A^ — 1 coupled equations at any order j € N. The first two orders (j = 0, 1) 
explicitly read: 



.(0) 



k-1 



(.,-..,)4°^-^E4V^ 



(41) 



<?=! 



■,(1) 



-i{ujk - kuji)v'j}'> + iA;i/(^)4°^ + 



.Wfc 



fc-i 



9=1 



(0) 



(42) 



Again, a sum is absent if the upper limit of summation is less than one. The latter two 
systems, in which the unknowns are vj,,vj, and u^^' , have a triangular structure, i.e. the 
equations can be solved one after the other starting from the first down to the last one, with 
the initial condition (j40p . The expansion of the frequency shift (|39p (i.e. the coefficients 
u''^') is determined by requesting the vanishing of possible secular terms in the expansion 
of vi. One thus arrives at the following leading order expansions: 



ih{fl) = UJl + 777^^ 7 + 



8(2^1 



^2 



(43) 



vi{t) = 1 



8(2a;i - ^2)^ 



^i{2uJi-uj2)t I _|_ 



(44) 
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V2{t) 



flU}2 



4{2uJi - U2) 



J{2ljJl—UJ2)t 



+ 



v?,{t) = — 7 I V 



8(a;2 — 2a;i) I oo^ — ijJ2 — o^i 

the dots standing for terms of higher order. 



uj^ — Soji 



+ ... 



(45) 
(46) 




Figure 11: e^ vs. t (log-log scale), k = 1,2,3,4, for FPU, E 
data. Black curves: analytic prediction. 



10" 



Red curves: numerical 



In Fig. [TT]the first four modal energies efc(t) of the a-model are plotted vs. time t, in 
log- log scale, at the energy value E = 10~^. The red curves represent the numerical values, 
while the four black lines correspond to the three analytic formulae (I44p -(H6]1 and to that 
of mode A; = 4, which is not explicitly reported here, since it has a very long, uninteresting 
expression. Note that e^ = \vk\'^. Note how the agreement is excellent up to the final 
computation time t = 10^. 

Fig. [TT] shows that the perturbative scheme used here works very well at low energies. 
The regime of validity of the perturbative solution can be estimated by requiring that in 
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the expansions (j43p and ()44p of f^i and of the vi, respectively, the leading order is dominant 
with respect to the correction proportional to fi'^. The latter condition turns out to be more 
restrictive than the former and, by taking into account that 2uji — uj2 = 7r^/(4A^^) + 0(A^~^), 
yields 



vr^ 



a'^EN^ < — ~ 10 , (47) 

8 

i.e. an upper bound to the energy which decreases with increasing A^. For A^ = 32 and 

a = 0.33, one gets E < 3 ■ 10~^. Using (|43p . condition (147p can be transformed into 



^"^ < ^ (48) 



2oJi — U2 2 

Thus the ratio of the nonlinear frequency shift 6uJi = Qi—loi to the resonance width 2uji—U2 
must be smaller than a constant, which is akin to the so-called resonance overlap criterion 
first formulated in Ref. [49] (notice that in the latter paper, dealing with the FPU-/3 model, 
the considered resonance width was different; see also Refs. [9] and [11] for comments on 
this point). 

For values of the energy larger than about 10^"^ the agreement of the expressions (|44p - 
(j46p with the numerical data starts to worsen, and one would really need to include higher 
order corrections. However, some leading order analytic predictions concerning the behavior 
of the VkS can still be obtained, within the same perturbative scheme, either on very short 
times, or at saturation, as follows. 

4.4 Secular avalanche 

Since by the initial condition (|40p all the vj,, are zero at t = if /c > 2 and ujk — kuii is small 
for low modes (and exactly vanishes for A; = 1), the contribution to the short term dynamics 
coming from the linear term (cjfc — kui)vj, on the right hand side of the equations ()4ip is 
expected to be small on short time. By neglecting such a term the leading order equations 
()4ip simplify to 

9=1 

where again the sum on the right hand side is absent for k = 1. Equations (j49p with the 
initial condition vj. = 6k^i admit the exact (unique) solution 

4°^ (t) = {-itf-^Ck (A; = 1, . . . , iV - 1) , (50) 

where the coefficients Ck are computed by solving the nonlinear map 

for 2 < fc < A^ — 1. Notice that to leading order, the first mode does not change from 
the initial condition: v\ ' {t) = 1. The evolution equation for the first mode at second 
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order can be obtained by solving the second order equations (j42p for k = 1 with the choice 
u''^' = 0.[50] Taking into account that v\ =1, and making use of (|50p and (j5ip to get 
V2 = —itc2 and C2 = (^2/4, one gets a differential equation for v\ , which can be integrated 
with the initial condition v\ = 0. The result is 

„«(*) = -^ . (52) 

Finally, the normalized modal energies e^ = |ffcp can be computed: 

e,(t) = |.inAr) + ...p:.l-^^^i^; (53) 



e.(0 = l/-^4°^ + ...P = (Mt)^('=-^)ci {k>2), (54) 

which explains the derivation of the law (|lip -(jl2 p reported in Section II and tested in Fig. [2] 
of Section III. 

4.5 Role of resonance detuning and energy conservation 

Let us discuss the validity of the approximation obtained by neglecting the linear term on 
the right hand side of equations (j4ip . This is done by substituting on the right hand side of 
the same equation the solution (j50p and by checking whether and how long the linear term 
stays smaller than the quadratic term. One thus finds that a first necessary condition for 
the approximation to be valid is that 

t^Tdet^k-J_ ^^^^2). (55) 

kuJi — Uk 

The detuning time T^'^^ is the characteristic time within which mode k is resonantly pumped, 
the term proportional to the resonance amplitude koji — uj^ being negligible. Now, T^"^* is a 
monotonically decreasing function of A;, and for high mode numbers T^^* ~ N , whereas for 
low modes T^^^ ~ N^ /k"^ , which can be as large as 0{N^) for the first few modes, actually 
too large to be compatible with the observed time-scale of formation of the quasi-state. 
When checking the validity of the approximation made by neglecting the linear term in 
()4ip . one has to take into account that the quantity /ii^^, defined in ()36p . is constant and 
equal to one. The latter condition implies that I'WfcP = fife < 1 for any k > 2, yielding, when 
the (approximate) solution (f50]l is taken into account, a second necessary condition for the 
validity of the approximation made, namely 

1 
r ''-'' 
t < T^" = -^ . (56) 

T^° is the mode dependent characteristic time of the secular avalanche process, i.e. the time 
within which mode k can grow, when resonantly pumped by the cascade process, without 
violating the conservation of the energy hi^-^ . We find that T^° is a monotonically decreasing 
function of k and, by direct inspection of the map (j5ip . that T|" ~ 0.2N/fi, if k is not too 
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small. Note that from §3^ it follows that ei > if t < Tf" ^ 2N/{-k^), while there is 
obviously no detuning time for the first mode. 

Having two characteristic times for each mode k > 2, namely T^"^* and T|'^, the actual 
characteristic time T^ at which e^ stops its secular growth is defined as the minimum of 
the two for each k: T^ = in.m{T^"',T^'^^}. For small values of /i, T^ = T™ for k < k^ and 
Tfc = T^"^* for k > kc- An order of magnitude estimate of kc is obtained by observing that 
T^" w N/ ji for almost all modes, whereas T^^* ~ N^/k'^ if k/N <^ 1 (numerical factors are 
neglected). Then the crossover r|" = Tj:'^^ takes place at the critical mode number 

^^VJ^, (57) 

valid if kc/N <^ 1, i.e. if // <^ 1. The time-scale for the formation of the quasi-state is 
defined as the crossover time at kc- 

which is the estimate (jlSp . 

Prom the above reasonings, it follows that modes with k < k^ are always inside a 
resonance layer and their modal energies stop to grow due to energy conservation, whereas 
modal energies e^ with k > k^. stop to grow due to resonance detuning. We thus distinguish 
core modes with k < kc from those with k > kc- Core modes share most of the energy in 
the quasi-state, and the resonant transfer of energy among them is effective. Higher modes 
outside the core are resonantly pumped for a certain time and are tuned out of resonance 
rather quickly. The quasi-state is characterized by a few core modes sharing the energy 
and by an exponential localization of the m.e.s. for high modes. Notice that for A^ = 32 
and a = 0.33, equation ()57p gives kc varying from 2 to 8 as the energy varies from 0.01 to 
1, in agreement with the number of modes that define the plateau of partial equipartition 
observable in the numerical spectra reported in the left column of Fig. [3l 

4.6 Fixed point of the averaged system: quasi-state 



Let us consider the averaged equations (p4|) . In the limit /i — )■ one has T-Ls — )■ 0, we require 
ill — ^ i^ii and the equations simplify to 

Vk = -i{ujk - kuji)vk , (59) 

whose solution with initial condition (|35p is 'yfc(t) = 'yfc(O) = <Ja:,i- In other words, the system 
is placed at t = on a fixed point and stays there forever. Notice that system (j59p admits 
a one-parameter family of fixed points 

Co = {\vi\ = 1 ; |t;fc|=0, k>2} , (60) 

which is a unit circle in the wi-plane of C^~^. It can be easily checked that the set Cq is the 
maximum of the Hamiltonian T-L2 — h^i (limit of the Hamiltonian ^ as // — )• 0) constrained 
to the set h^)-^ = 1. If /U > but small enough, the set Cq is slightly deformed to the set 
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Cfj_ of stationary points of the full averaged system (j34p on the set h^-^ = 1, i.e. the set 
of critical points of the Hamiltonian A = ^2 — /if^i + ^3 constrained to the set h^^ = 1. 
Moreover, if /i is small enough, the critical set C^ preserves its maximum property (i.e. C^ 
is the maximum of A on h^^ = 1) and turns out to be Lyapunov stable: if v is initially 
close to C^ it stays close to it forever. On the other hand, the initial condition ffc(O) = 5k,i 
means that one starts on Co, so that if Co and C^ are close, i.e. if /i is small enough, the 
orbit of the averaged system stays forever close to C^. We make use of this qualitative (but 
exact) reasoning to identify the quasi-state orbit of the a- model with the set C^. 

The fixed points of system (I34p can be perturbatively computed by means of the ex- 
pansions (J38|) - (j39|) . To leading order, one has to find the fixed point of system (j^T]) . setting 
VI , which yields 



fc-i 
(0) , (0) Wfe >r^ (0) (0) 



^ vrA> ' (^ ^ 2) (61) 



where A, at this stage, is an arbitrary complex number such that |Ap = ei, to leading order. 
Isolating the dependence of t;^ on A one gets 

vf^ = X'fk , (62) 

where the real coefficients fk are found by solving the nonlinear map problem 

fc-i 

/l = 1 ; fk = 777 7 Y] fk~qfq (63) 

A{kui - Wfc) ^ 

for 2 < A; < A^ — 1. Recalling that, to leading order, v^ = pf^~^v), , one can compute the 
m.e.s. 

efe = \vk? = fl (/^^ei) " ei . (64) 

Such a spectrum depends on the unknown ei = |Ap, which is determined by the normal- 
ization condition h^^^ = 1, i.e. 

'.«^E(^)-^E(^)/i(.%re.^- (») 

The latter condition determines the modulus of A, so that an overall phase remains unde- 
termined. For this reason, the set C^ is composed by a one-parameter family of stationary 
points of the averaged system. 

We now observe that to a stationary solution of the averaged system, there corresponds a 
periodic orbit, or one-dimensional torus, of the original system: recall that, before averaging, 
we performed the change of variables u^ = e~^^^^^Vk- The fundamental frequency ili of the 
quasi-state periodic orbit is easily determined by means of equations (|12|) . Setting there 
k = 1 and u} = 0, and making use of ([62]) - ([63]) . one gets 

(0)/ (0)x* 

2yf) 8(2^1 - ^2) ^ ^ 
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which imphes 

n , 2 (1) , /i^eia;iL^2 .^„s 

ill = ui + fi-'u^'' = bJi + —- r- (67) 

8(2a;i - uj2) 

Notice that in the regime of strong locahzation, when ei ~ 1, the above expression for the 
corrected frequency coincides with that given in (j43p . 

As already mentioned, the quasi-state m.e.s. (j64p fits the numericahy obtained one 
at fow energy, when the exponential localization is so strong that only a few low modes 
contribute to the total energy. In this case, the map ([63j) can be approximately solved, as 
follows. First of all, notice that if only low (acoustic) modes are relevant, one can substitute 
the approximate expressions uj^ ~ Trk/N, kuji —cok — iT^k{k'^ — 1)/(24A^^) in the numerator 
and denominator of (j63p . respectively. Second, we look for a solution of the map in the form 
fk = k g , which fits the condition /i = 1. Upon substitution of the latter expression into 
()63p . and after some simple calculation, one gets g = vr^/A^^. Thus, in the regime of strong 
localization, the approximate solution of the map (j63p is given by fk — k{N'^ / tt'^)^~^ . As 
a consequence, recalling that /i^ = o?E jN and that ei = Ei/E, for the spectrum ()64p one 
gets the approximate expression 

^k = k'{ :^\ ei . (68) 

Under the same acoustic approximation (w^ ~ kuJi) leading to the latter formula, condition 
([65]) becomes 

N-l N-l , 2 7^Ar3 \ fc-1 

E'-E'^^^ , = 1. ,69) 

k=\ k=\ ^ ^ 

This explains formulas (J16p -(|18p given in Section II. We stress again that the m.e.s. (|68p 
corresponds to the leading order spectrum of the Lyapunov continuation to the nonlinear 
case of the first mode, first reported in Ref. [5] and there named g-breather. We also 
recall here that the normalization condition (j69p can be simply solved. Indeed, by defining 
X = a^EN^ /tt'^ and y = a^ EiN'^ / tt'^ , it is easily found that (|69p . in the limit of large A^ and 
for y < 1, yields the simple implicit equation y(l + y)/{l — y)^ = x. This is explicitly solved, 
for X > (and < y < 1), to give y = F{x), for a suitable function F, which determines 
El (and thus ei) in terms of the total energy E, of a and of the number of particles A^. In 
particular, having determined F{x), it turns out that, at any fixed x, ei = F{x)/x, which 
can then be used to finally determine the m.e.s. (|68p . 

To conclude this subsection, we observe that all the above construction was performed 
to leading order, and applied to the Toda model, would give the same result. This is a 
consequence of the coincidence to third order of the FPU-a and the Toda models expressed 
by relation @. 

4.7 Long term stability of the quasi-state 

The argument used to guess that the one parameter family C^ of stationary points of the 
Hamiltonian (j33p constrained to the set h^j^ = 1 is the candidate quasi-state orbit of the 
FPU-a and Toda models, is a qualitative one. We did not specify how small the parameter 
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/i has to be to ensure that C^ is a inaximum of the constrained averaged Hamiltonian and 
moreover the unperturbed orbit Cq (which contains the initial datum) is in the domain of 
concavity of C^. In this respect, we do not provide any precise estimate on fi, which would 
require technical arguments that are beyond the aim of the present paper. We only mention 
that, from the discussion made in Section IV C above, one may expect that a condition of 
the form E = 0{N~'^) should hold to ensure the Lyapunov stability of C^. 

Apart from the precise energy-size dependence of the stability domain of the one- 
dimensional torus C^, it is important to understand that the construction made above 
is a leading order one, which ensures, for the real system, confinement close to the set C^ 
on a relatively short time-scale, typically an inverse power law of the small parameter ^. 
In order to explain the stability on much longer times observed numerically, one should go 
beyond the leading order construction, which modifies the averaged Hamiltonian (|33p with 
the addition of higher order terms, and increases the length of the confinement time. Such 
an approach was used in Refs. |5H [52} [53] to prove that in the case of initial excitation 
of the first mode in a suitable class of string equations, the confinement time grows as the 
exponential of the inverse of the small parameter. However, it must be stressed that im- 
porting those techniques to the problem at hand is nontrivial at all, essentially due to the 
lack of complete resonance in the FPU problem (explicit reference to the FPU problem, in 
this perspective, has been made in Ref. |52j). 



5 Concluding remarks 

We have shown that the classical FPU paradox can be easily understood if one compares 
the behavior of the (presumably) nonintegrable FPU-a-model with that of the integrable 
Toda model, for the same initial condition, at the same energy. One thus finds that the 
dynamics of the a-model on a short initial time interval is very close to that of the Toda 
model. The quasi-state guessed by FPU is nothing but the quasi-periodic Toda torus on 
which the dynamics is initiated for the same initial condition. At very low energy, such a 
torus becomes one-dimensional and the FPU-a system evolves close to it up to extremely 
long times. As the energy increases, the dimension of the reference torus grows, while the 
time of confinement close to it (or stability time) for the FPU-a system decreases. Finally, 
the approach to a state of energy equipartition of the FPU-a model takes place through a 
diffusive-like detachment from the quasi-state, a phenomenon that is chaotic in character, 
and consequently, absent in the integrable model. 

Regarding the FPU a-model as a perturbation of the Toda one, as first suggested in 
Ref. |16j . and strongly motivated here, seems to be quite promising. For example, in Ref. 
|54| it is shown that the Toda Hamiltonian fits the hypotheses requested to the unperturbed 
Hamiltonian in the Nekhoroshev theorem, with the purpose to explain the phenomenology 
of the FPU-a model in terms of the exponential stability of all the Toda integrals (or 
actions). However, precise/optimal estimates of the regime of validity (energy-size domain) 
of such a conclusion are still lacking. Of course, if the energy is small enough, one cannot 
exclude the possibility to enter a KAM regime, i.e. to have stability over infinitely long 
times. Here also, we know that this is possible in principle [55], but again, precise/optimal 
estimates of the regime of validity are not available yet. We finally notice that it is not 
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presently clear whether this approach (i.e. considering the FPU system as a perturbation 
of the Toda one) is valid for initial conditions involving also high frequency modes or even 
for generic excitations. 

One interesting question is whether and how the energy region for the long-term stability 
of the quasi-state collapses or not to zero by increasing N. In the light of the purely 
numerical results contained in Ref. [28], a first answer seems to be that exponentially 
long times to equipartition might be, for the a-model, a finite-size effect, observable only 
in an unspecified regime E ^ as N ^ oo. In the same paper, it is shown that for 
particular low frequency initial data, the /3-model seems to exhibit exponentially long times 
to equipartition in the thermodynamic limit. On the other hand, in Ref. [56], an analytic 
lower bound to the canonical mixing time of a weakly coupled (j)"^ model is provided, showing 
that such time is longer than a stretched exponential of the inverse of the temperature. The 
full scenario is quite complicated and deserves a lot of further study. 
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